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We evaluate the degree of quantum correlation between two fermions (bosons) subject to contin- 
uous time quantum walks in a one-dimensional ring lattice with periodic boundary conditions. In 
our approach, no particle-particle interaction is considered. We show that the interference effects 
due to exchange symmetry can result into the appearance of non-classical correlations. The role 
played onto the appearance of quantum correlations by the quantum statistics of the particles, the 
boundary conditions, and the partition of the system is widely investigated. Quantum correlations 
also been investigated in a model mimicking the ballistic evolution of two indistinguishable particles 
in a ID continuous space structure. Our results are consistent with recent quantum optics and 
electron quantum optics experiments where the showing up of two-particle non-classical correlations 
has been observed even in the absence of mutual interaction between the particles. 

PACS numbers; 05.40.-a, 03.65.Ud, 03.67.-a 
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I. INTRODUCTION 

Quantum walks (QW) describe the random walk be- 
havior of a quantum particle Due to quantum me- 
chanics effects, such as the coherent superposition of 
wavefunctions and interference, QW exhibit a qualita- 
tively different behavior with respect to classical random 
walks, such as the ballistic propagation of the wavefunc- 
tion instead of the diffusive behavior exhibited by the 
classical probability distribution [2|]. Quantum effects 
make QW extremely promising for the implementation 
of more efficient and faster research algorithms than the 
protocols commonly adopted in the classical computa- 
tion [^-Q . Due to their potential application in quantum 
information science, simple models of QW have been in- 
vestigated and physically implemented in various physi- 
cal systems ranging from quantum optics ^6j-£8|] to nuclear 
magnetic resonance setups d, • 

Two kinds of QW are considered in the literature: 
discrete-time quantum coined walks and continuous-time 
quantum walks (CTQW). In the former, a two- level state, 
the so-called coin, rules the unitary discrete-time evo- 
lution of a particle moving in a lattice of sites. On 
the other hand, for CTQW the evolution of the parti- 
cle is continuous in time and it is only determined by 
a Hamiltonian whose terms represent transitions among 
the lattice sites. CTQW have successfully been imple- 
mented in lattice waveguide systems jOj-Q , where the ap- 
pearance of non-classical correlations has been observed 
for two-photon input states [1, Indeed, due to in- 
distinguishability of the photons, probability amplitudes 
of two-particle wavefunctions can interfere, thus leading 
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to the formation of correlations with no particle-particle 
interaction. Interference effects have also been success- 
fully detected in the so-called electron quantum optics ex- 
periments, namely quantum optics experiments-like with 
electrons in solid state systems, where many-particle ef- 
fects play a key role [HI, Recent developments in 
nanofabrication technology have allowed exploiting the 
ballistic electron transport along chiral edge states at 
integer quantum Hall regime for the experimental re- 
alization of Mach-Zehnder and Hanbury Brown- Twiss 
interferometers [ill, [13 • these studies, two-particle 
Aharanov-Bohm oscillations, due to exchange symmetry, 
have been related to the degree of quantum correlation 
between the charge carriers, even if electron-electron in- 
teraction is not taken into account [Tsj . 

In the last decade, the notion of quantum correlations 
(QC) as a fundamental source for the implementation of 
quantum information algorithm and commonly referred 
as entanglement has been widely investigated in systems 
of identical particles [T^ - flsj . Indeed, its quantification is 
certainly crucial for understanding a number of physical 
phenomena involving correlated indistinguishable sub- 
systems. The main difficulties appearing in the definition 
of a criterion able to classify and quantify the amount of 
QC among indistinguishable particles are closely related 
to exchange simmetry which requires the symmetriza- 
tion or the antisymmetrization of the quantum wavefunc- 
tions describing bosons or fermions, respectively. Differ- 
ent approaches have been used to estimate the degree 
of non-classical correlation in bipartite systems of identi- 
cal particles [T^ - [l6j . The Schliemann criterion relies on 
the fermionic analogous of the Schmidt decomposition, 
namely the Slater decomposition pjjj. In the approach 
developed by Zanardi [15], the entanglement is evaluated 
in terms of the QC between modes by mapping the Fock 
space of the modes themselves into qubit states. In the 
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criterion proposed by Wiseman and Vaccaro [T^, ITsj , the 
entanglement of the particles is a sort of accessible en- 
tanglement, i.e. the maximum amount of non-classical 
correlations which can be extracted from the system by 
means of local operations and then transferred into con- 
ventional quantum registers. 

While the correlation between position and coin de- 
grees of freedom of particles has been widely investi- 
gated in discrete-time quantum walks [H, Hs-HH, an 
exhaustive analysis of the building up of the QC in 
CTQW of identical particles is still lacking. The appear- 
ance of non-classical correlations in two-particle quan- 
tum optics CTQW setups has been discussed only qual- 
itatively 6, 8j. The aim of this work is to provide an 
accurate analysis of the quantum correlation created in 
the CTQWs of two identical particles (both fermions 
and bosons) in a one-dimensional (ID) system. Specifi- 
cally, we will first examine the diffusion of two fermions 
(bosons) in a lattice of sites with periodic boundary 
conditions, in order to study the development of non- 
classical correlations in simple models mimicking exper- 
imental quantum optics and electron quantum optics se- 
tups. Then, we will analyze such a system in the limit of 
a large number of lattice sites, with a vanishing intersite 
distance, thus resulting into the free propagation of two 
identical particles along a ID structure. 

Our model does not take into account any kind of 
particle-particle interaction, so that the time evolution of 
the two-particle quantum state involving exchange sym- 
metry is essentially ruled by single-particle Hamiltonians 
describing CTQW. In this way, the appearance of QC 
is only related to two-fermion (-boson) amplitudes in- 
terference due to the quantum statistics of the particles 
involved in the process. In this work, in order to estimate 
the degree of non-classical correlation between the posi- 
tions of the particles, we adopt the criterion proposed by 
Wiseman and Vaccaro In fact, such a criterion, 

unlike the approach developed by Schliemann fl^, be- 
haves correctly under one-site (local) and two-site (non 
local) transformations. Furthermore, its use does not 
lead to the violation of the local number of particle super- 
selection rule, as may happen in the case of the procedure 
proposed by Zanardi [15| . 

The paper is organized as follows. In Sec. |lTl we de- 
scribe the theoretical approach used to quantify the QC 
between two indistinguishable particles. In Sec. IIIIl we 
study the CTQW of two fermions (bosons) in a lattice 
of sites with periodic boundary conditions and then es- 
timate the time evolution of the amount of QC of the 
system. In Sec. lIVl we numerically evaluate the degree of 
non classical correlation between two identical particles 
free propagating in a ID structure. Finally, conclusions 
and discussions are given in Sec. |Vl 



II. EVALUATOR OF THE DEGREE OF 
QUANTUM CORRELATIONS IN BIPARTITE 
SYSTEMS OF IDENTICAL PARTICLES 

Here, we briefly illustrate the theorethical criterion 
adopted to evaluate the degree of quantum correlation 
in two-fermion and -boson systems. It is based on the 
notion of entanglement of particles proposed by Wise- 
man and Vaccaro [H, [ll] . 

By using the mode-occupation representation, an ar- 
bitrary pure two-fermion (boson) state in a M-mode sys- 
tem can be expressed as 

l*) = E^wlM), (1) 

where the integers rii of the set {n} = ni, . . . , n^, . . . , um 
satisfy the relation ni -1- . . . -f -I- . . . + um = 2. Here, 
the ket \{n}) indicates the state vector in the Fock space 
with TiiS particles in the i-th mode, and the C{„j.'s are the 
coefficients of the linear superposition. While for bosons 
rii's range from to 2, for fermions the occupations num- 
bers are restricted to be or 1 due to the Pauli exclusion 
principle. 

As argued by some authors ^IS^, a formal equivalence 
between the space of the occupation-number states and 
the tensor product space of the modes can be established. 
In this way, the occupation number of each mode consti- 
tutes a distinct state of the mode itself. In the Zanardi 
approach [l^, the amount of non-classical correlations 
between the occupation numbers of the modes controlled 
by two parties of the system, namely Alice and Bob, rep- 
resents the so-called entanglement of modes. The latter 
does not always constitute a valid measure of the true de- 
gree of quantum correlation between Alice and Bob [l^ ■ 
As a matter of fact, not only the entanglement of modes 
can give values different from zero even when applied to 
suitable single-particle states [l^], but also it does not 
take into account the local-particle number superselec- 
tion rule (LPNSR)[16]. Indeed to fully exploit the QC 
between modes, Alice and Bob, at least in principle, must 
be able to arbitrarily measure and manipulate their local 
systems. Unless each party of the systems possesses a 
definite number of particles, this will lead to a violation 
of the LPNSR. 

The Wiseman and Vaccaro criterion satisfies the LP- 
NSR. In such an approach, in addition to the two iden- 
tical particles shared by Alice and Bob, their quantum 
state 1^) involve a standard quantum register, namely a 
set of distinguishable qubits. The entanglement of the 
particles Ep \s defined as the maximum amount of QC 
that Alice and Bob can produce between their standard 
quantum registers by means of local operations. As a 
consequence of the LPNSR, i?p, in place of the QC be- 
tween the modes that Alice and Bob have access to, is 
given by 

2 

£;p(|M')) = E/'„i?M(|*("^)), (2) 

n=0 
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where |\I>(")) is the projection of the quantum state j^*), 
describing the global system, onto the Fock subspace 
where Alice controls n particles and Bob the remaining 
(2-n) ones. P„ = (^'(") is the probability for Alice 

(Bob) of finding n {2 — n) particles as a consequence of 
a measure of the local number of particles, while Em is 
the degree of quantum correlation between the two sets 
of modes, each controlled by a party of the system. In 
other terms, Ep is the weighted sum of the entanglement 
of modes when the local particle number is measured. It 
is worth noting that Ep depends upon the partition of 
the system, that is upon which modes Alice and Bob con- 
trol. From this point of view, different partitions of the 
system can lead to different values of Ep. 

For two-particle pure states, the expression given in 
Eq. ^ takes a simple form. Indeed, only some quantum 
states 1 5'^"'), belonging to the Fock subspace with a fixed 
local number of particles, give a non vanishing contribu- 
tion to entanglement. Specifically, let us consider the 
case of the set of states l^f^"^) where the local number of 
particles possessed by Alice is zero. Any state of such a 
set can be written, in the mode occupation number, as 
|0„^) (S> |2„^) and it is separable, that is it can be factor- 
ized in a term describing Alice with no particle and in a 
term describing Bob with two particles. This implies that 
the contribution of to Ep is zero. Analogously, 

the amount of quantum correlations stemming from set 
of states I^E'^^)) is also zero. This implies that the bipar- 
tite entanglement is non-vanishing only when both Alice 
and Bob have one particle, that is only for the quantum 
states belonging to the set j^'^^^). As shown in Ref. p^ . 
Ep becomes 



(3) 

where e indicates an arbitrary quantum binary entropy 
and 
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is the single-particle matrix describing the subsystem 
controlled by Alice. The latter is obtained from l^'^^^) 
by means of the creation (annihilation) operators aj^iak) 
acting on Alice's modes only. The form of Ep given in 
Eq. ([3]) holds for either bosons or fermions. 



III. TWO-PARTICLE CONTINUOUS-TIME 
QUANTUM WALKS ON ID RING LATTICES 

In this section, a model of two-fermion and two-boson 
CTQW in ID graphs is solved by means of analytical 
and numerical techniques for a small and large number 
of nodes, respectively. The time evolution of the wave- 
function describing the system is then used to estimate 
the degree of quantum correlation according to the crite- 
rion given in Eq. Q. 

Here we examine CTQW on a ID ring lattices of 
N sites (with N even) with periodic boundary condi- 



tions. I n ag reement with previous single-particle investi- 
gations [23, the topology of the graphs here considered 
is simple in the sense that each node is connected to 
its two first neighbours. Even if more complex networks, 
among which two-dimensional lattices [IJ] or graphs with 
larger connectivity [25| . could be examined, our model is 
good enough to describe experimental implementations 
of CTQW, such as two-photon transport in an array of 
waveguide lattices where non-classical correlations ap- 
pear. Therefore, though simple, it represents a valid tool 
to analyze the amount of QC appearing in two-particle 
CTQW. 

The two-particle Hamiltonian describing the dynami- 
cal evolution of the system is given by 
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n — rt^ -f- ftp, 



(5) 



where ^^(^) is the single-particle Hamiltonian acting on 
the particle a{f3): 

'^"(/3)li>a(« = 7 (2b)Q(/3) - \ j - l)a(/3) - \ j + l)a(/3)) ■ 

. . (6) 

\j)a{fi) indicates the quantum state describing the par- 
ticle a(/3) localized in the j node and forming a com- 
plete, orthonormalized basis set, which span the whole 
accessible Hilbert space. 7 denotes the intersite trans- 
mission rate. Due to the periodic boundary conditions, 
here we assume that the node iV -I- 1 coincides with node 
1. Eq. ^ is the discrete version of the Laplacian, and, in 
turns, the discrete version of the Hamiltonian describing 
the free propagation of a particle in a lattice. 

Given the form of the two-particle Hamiltonian in 
Eq. ([5]), single-particle dynamics allows one to estimate 
the time evolution of the two-fermion and -boson system 
in the basis states given by 



(7) 



and 



I \j3)ap ■■ j = k, 

respectively. Specifically in order to solve the model, 
which can be interpreted as the free propagation of two 
identical particles in a periodic system, we adopt the 
Bloch function approach, in analogy with the method 
commonly used in solid-state physics [2^. The single- 
particle Bloch states 



(9) 



are eigenstates of the Hamiltonian of the Eq. © with 
eigenvalues En = 27 (l — cos ^^), and can be used to 
evaluate the coefficients describing the transition ampli- 
tude \k,jit) = {k\e-'"°*\j) from the state \j) at t = to 
state \k) at time t. Indeed the latter can be expressed as: 
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exp {—2ijt) 
iV 



E/ / 2n7r\ \ / 2Trn 

exp [ 2i7i cos [ — I j exp [-1 — (fc-j) 



(10) 



where is set H = 1. 



By using the two-particle basis states written in 
Eqs. and ([7]), we can evaluate, in terms of the above 
single-particle coefBcients, the two-boson (-fermion) 



transition amplitude fi^^''\g .jj.{t), from the state \jr)f(b) 
at i = to state \ks) /(t) at time t. It reads: 

f^^ks,jr{t) = Afc,,(t)A,,,(i) - \kAt)Kjit), (11) 

and 



Xk,jit)Xs,rit) + Xk,r{t)Xs,j{t) k^s and j^r 
/fcsjrW = <! V2Afc,,(t)A«,,(t) fc = s xor j = r . (12) 

Afej"(i)As.r(0 k — s and j = r 



Once the time evolution of the two-particle quantum 
state \jr)f{b) describing two fermions (bosons) initially 
localized in the node j and r is known, the dynamics of 
the built up of the amount of quantum correlation can be 
quantified by using Eq. ([3]). To this aim, different parti- 
tions of the system can be considered. Note that here the 
modes of the system correspond to the sites of the lattice. 
In the simplest case, Alice controls the first half of the lat- 
tice and Bob the second one, namely A = {1, . . . , -^l and 
B = -f 1, . . . , N}. Thus, the single-particle density 
matrix of the Alice's subsystem given in Eq. Q takes, 
for the initial state the form 

N 



E A*Ljr Wife's jr(i) 



kk' ^ N 



(13) 



Its quantum entropy can be quantified by means of the 
von Neumann entropy 

e„^(t) = -Tr{[p(t)W]ln[p(t)^^)]}, (14) 

which represents an appropriate evaluator of the degree 
of the correlation in standard bipartite systems, and this 
can be used in Eq. ([3]) in order to estimate of the entan- 
glement of particles. For sake of completeness, in the fol- 
lowing, other possible partitions of the system will be ex- 
amined, where Alice and Bob can access to non-adjacent 
modes. The evaluation of Ep in such cases requires to 
calculate again both the single-particle density matrix of 
the subsystems and the probability of finding one parti- 
cle in each subsystem, being the latter closely dependent 
upon the specific partition. 



A. Small number of nodes 

Here we study, by means of analytical techniques, the 
quantum correlation dynamics for the case of an ID ring 
lattice composed of a small number of sites. Specifically, 
we take iV = 4. Various initial configurations of the 
bosonic (fermionic) system, that is the states describing 
at time t=0 the two identical particles, should be exam- 
ined. However, due to the periodic boundary conditions 
of the lattice under investigation, the number of initial se- 
tups of the system resulting into different time evolutions 
of the QC turns out to be smaller than the dimension of 
the Pock space of the two-fermion (-boson) system. 

First, let us analyze the case of particles initially local- 
ized in two different nodes for the partition A — {1,2} 
and B — {3,4}. Three configurations are possible: both 
particles occupy the Alice (Bob) modes, that is the quan- 
tum state |12)j:(f,) (|34)j^(ft)); one particle is in an Alice 
site and the other one is localized in the adjacent site con- 
trolled by Bob, the state |23) or 1 14) y (t) , or in the non- 
adjacent Bob mode, the state |13)y(h) or |24)y(j). Thus, 
to analyze the role of the initial configuration of the sys- 
tem into the building up of non-classical correlations, we 
just need to evaluate the time evolution of the quantum 
states |12)y(b), |23)/(fc), and |13)/(;,). As shown in the 
previous section, the latter can be calculated from the 
single-particle amplitude transitions given in Eq. (jlOl) . 
In agreement with the results of Ref. [23], these, after a 
straightforward calculation for the case of A^=4, can be 
expressed as: 

Xk,k {t) = exp (-2z7t) cos^ jt 
Xk±i,kit) = Afe,fe±i(t) = exp {-2ijt)i sin "ft cos "ft 
Xk±2,k (t) = Xk,k±2 (t) = - exp {-2ijt) sin^ jt. (15) 
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By inserting the above coefBcients in the Eqs. p2|) and we can quantify the building up of QC for the initial 
PT|) . and then by using the expressions ([T^ and ([T^. quantum states of interest. We find that 



fib) 



A={1,2} 



[cos(27 i) - 1]2 , r [cos(27 t) - if 
: In 



[cos(27 t) + iX" 



In 



2[l + cos2(27 t) 
[cos(27 t) + 1]2 1 



2[l + cos2(27 t)] 



(16) 



In Figure [U we report the time evolution oi Ep, namely 
the amount the QC, for a two-fermion (-boson) system 
prepared in different initial conditions. For the sake of 
simplicity, in this section we call time the adimensional 
parameter jt. In agreement with theoretical predictions, 
at t—0, Ep is null, thus indicating the absence of QC 
between the two particles localized in different nodes. 
After initial time, Ep remains zero for the input state 
\^^^^^) fib)^ while it exhibits oscillations, between and 
In 2/2, with a period of 7r/2 for the other two input states 
/(b) and \ip^^^^) jj-b) . While in the former initial con- 
figuration the dynamics of the two-particle wavefunction 
does not result into building up of non-classical correla- 
tions, in the latter cases the interference effects are able 
to produce QC, whose periodic behavior is strictly re- 
lated to the evolution of the system in a ring lattice with 
a small number of nodes with periodic boundary condi- 
tions. Here, the highest value of Ep is In 2/2, that is 
the half of the maximum degree of quantum correlation 
achievable between two particles, each one in a two-mode 
subsystem. It is worth noting that, for the partition con- 




FIG. 1. (Color online) Time evolution of Ep for three input 
states, namely |12)j^(t.) (solid line), |13)/(f,) (dashed line), with 
Alice controlling sites {1, 2} and Bob {3, 4}. For sake of clar- 
ity, the entanglement dynamics of the initial state |23) j^(t.) has 
not been reported being always null. 



sidered, Ep turns out not to depend upon the quantum 
statistics of the particles involved in the process. Such a 
result is somehow unexpected due to the different num- 
ber of states accessible to the two systems during the 
time evolution. Indeed, unlike fermions, bosons can oc- 
cupy simultaneously the same node, that is the transition 
probability to the states \kk)bS does not vanish as time 
goes by. Even if the time evolution of the bosonic wave- 
function differs from the fermionic one, for the specific 
partition here considered the two-particle interference ef- 
fects lead to the same single-particle features, as, in par- 
ticular, the reduced density matrices of the subsystems. 
As a consequence, also the amount of QC is the same. 

The results above reported are closely related to how 
the nodes of the ring lattice are assigned to Alice and 
Bob. If we move from the partition A = {1,2} and B = 
{3,4} to the one A = {1,3} and B = {2,4}, not only 
the QC turn out to depend upon the quantum statistics 
of the particles, but the highest degree of correlation can 
also be reached. Let us examine the time evolution of two 
bosons, or two fermions, initially occupying the nodes 1 
and 3. After a straightforward calculation, we obtain 
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Epdip'-'^^hf] =ln2sin2(27t) and Ep (iTp'-^^hb] =0. (17) 

V y A={1,3} V /A={1,3} 



While in the two-boson system no correlation is cre- 
ated, in the two-ferniion one Ep exhibits periodic oscil- 
lations and reaches the maximum value, namely In 2, for 
7* = {2k + l)/(47r) with keN. 

As stated above, Bose-Einstein statistics allows for the 
localization of two bosons on the same site. When the 
latter is taken as initial configuration of the two-particle 
system, the amount of non-classical correlations built up 
does not change with time and remains equal to zero for 
any partition examined. This means that time-evolved 
states of the kind \kk)b can always be factorized in terms 
of an one-particle state in the Alice modes and one in the 
Bob modes. Thus, an initial bosonic bunching prevents 
the appearance of QC in CTQW. 



B. Large number of nodes 

For a lattice with N sites, the dimensions of the Hilbert 
space for the two-boson and two-fermion systems are 
N{N + l)/2 and N{N ~ l)/2, respectively. This makes 
the analytical techniques inefficient for large N both to 
solve the two-particle dynamics and to evaluate the de- 
gree of quantum correlation. For this reason, we imple- 
mented a numerical approach that, once known the time 
evolution of the state, allows one first to diagonalize the 
single-particle reduced density matrix given in Eq. (|13p. 
and then to estimate, at any time, Ep. Such an approach 
is used here to quantify the amount of quantum correla- 
tions stemming from two-particle CTQW in an ID ring 
lattice with N — 50, 70, 100, when Alice accesses the first 
half of the lattice and Bob the second half (namely, the 
partition A = {!..., N/2},B = {N/2 + l,...,N}). 




FIG. 2. (Color online) Time evolution of Ep for the CTWQ 
of two fermions and two bosons, in the left and right panel, re- 
spectively, on a lattice with a number of sites N = 50 (dashed 
line). A'' = 70 (solid line), and iV = 100 (dotted line). Both 
systems are initially prepared in the state \ yy ~^ ^)f(b)- 

Fig. [2] displays the time evolution of Ep for two 
fermions (bosons) initially localized in two adjacent sites 
belonging to different subsystems, corresponding to the 
input state | A^/2, A^/2 + 1) As expected, at i = the 
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amount of QC is zero. As time goes by, Ep increases 
and saturates, apart from small oscillations, around a 
specific value for a time interval 7T which is found to be 
linearly dependent upon the number N of lattice nodes. 
These oscillations are closely related to the two-particle 
interference effects stemming from the spreading of the 
spatial wavepackets in a discrete lattice. In the limit of 
large values of N, the fluctuations disappear leading to 
a smoother behavior, as it will be shown in the next sec- 
tion. After a time 7T, Ep shows first a rapid increase 
and then again large fluctuations around a new satura- 
tion value. 

In order to get a better insight into the appearance of 
QC, we focus on the time-evolution of the two-fermion 
(-boson) correlation function ^j.^^'^t) namely the 

square modulus of the coefficient n .At)- To this 

aim, we report in Fig.[3]the values of F^j'''' (t) , at four dif- 
ferent times, for the case of a ring lattice with 70 nodes. 
At initial time, the two fermions (bosons) are localized 
in two adjacent sites (taken in the middle of the lattice) 
in a separable state. After a relatively short time, their 
wavepackets overlap and interfere thus leading to the first 
rise of Ep. Then the fermions (bosons) freely propa- 
gate along the lattices exhibiting spatial antibunching 
(bunching) S6f], as shown by the two peaks positioned 
along the antidiagonal (diagonal) of the correlation ma- 
trix (see panels (b) and (c) of the Fig. [3]). During this 
ballistic evolution, the amount of QC exhibits small fiuc- 
tuations around a given value in the interval 7T. The 
latter can now be viewed as a characteristic parameter 
of the dynamics of the system. Indeed, it corresponds to 
the travelling time of the two particles along half a lat- 
tice before the transmission thorugh periodic boundary 
conditions occurs. As expected, 7T is independent upon 
the quantum statistics of the particles. At longer times, 
the two-particle wavefunction is transmitted by the peri- 
odic boundary conditions, and their components interfere 
again creating additional QC. As reported in the panels 
(d) of Fig. |3]), at 7^ = 40, F^j''^ reveals a checkered pat- 
tern where both fermions and bosons are scattered along 
the whole lattice and Ep oscillates in time. 

The time evolution of Ep deserves further comments. 
The results reported here differ significantly from the 
ones found for the case of a lattice with = 4. First, 
QC depend now upon the quantum statistics of the the 
particles considered, though the qualitative behavior of 
Ep of the two-fermion and -boson systems is the same. 
Specifically, we find that the degree of correlation be- 
tween two bosons is always equal or greater than the one 
of the fermionic system. Such a behavior is at a first anal- 
ysis counterintuitive . Unlike the two-fermion dynamics, 
the bosonic bunching lowers the probability Pi of find- 
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ing each of the two particles in a different subsystem 
(see Fig. H]). As a consequence, according to Eq. ([3]), the 
amount of QC between two bosons should be lower than 
the one created in the fermionic system. However, also 
the degree of non-classical correlation of the two-particle 
state describing one particle in the Alice modes and the 
other one in the Bob modes, and quantified by the von 
Neumann entropy, has to be taken into account. In this 
case, the latter is larger for the two-boson CTQW and 
the product Pie„jV: namely Ep, turns out to be larger 
than the one corresponding to the two-fermion system. 

Furthermore, we find that the behavior in time of Ep 
is not periodic. Indeed, as expected, for lattices with 
a large number of sites, the transmission through the 
boundary conditions does not affect very much the time 
evolution of the two-particle wavefunction. As a conse- 
quence the interference among the different components 
of the wavefunction leads to a non-periodic behavior of 
the correlations Thus, the effect of the lattice periodic- 
ity on two-particle dynamics and the degree of quantum 
correlation is lower for lattices with higher values of N 
and should vanish in the limit of a very large N. 

IV. TWO-PARTICLE FREE PROPAGATION IN 
ID SPACE-CONTINUOUS STRUCTURES 

In this section, we numerically study the building up of 
non-classical correlations in the free-propagation of two 
identical particles in a ID structure. 

The physical model here investigated consists of two 
non-interacting massive fermions (bosons) travelling in a 
space-continuous structure. The latter can be viewed as 



with the normalization condition 

J^^ dxi J^^°° dx2\^b(f){xi,X2,to)\ = 1. The vari- 
ance a of the two single-particle Gaussian wavepackets 
and the distance |2xo| between their centers are 
such that, at the initial time, their spatial overlap is 
practically null. 

Here we estimate the amount of QC between the spa- 
tial degrees of freedom of the two particles. Specifically, 
we assume that the space domain [— cx), 0] is controlled 
by Alice, while Bob controls to the spatial modes in the 
interval [0,-|-oo]. To quantify Ep we adopt a numeri- 
cal procedure. In fact, even if the dynamics of the sys- 
tem could be evaluated analytically, nevertheless the es- 
timation of the time evolution of Ep as given in Eq. (|3]) 
would require involved calculations. On the contrary, 
the implementation of a numerical approach, exploiting 
the Crank-Nicholson finite difference scheme to solve the 
two-particle time-dependent Schrodinger equation, per- 
mits to evaluate efficiently, at fixed time steps, the built 



an infinite lattice with vanishing intersite distance with- 
out the interference effect stemming from the transmis- 
sion through the periodic boundary conditions. In other 
terms, we are considering the continuous-space limit of 
the model examined in the previous section. If the 
free propagation of massive particles obeying the Bose- 
Einstein statistics constitutes a helpful theoretical tool 
to investigate the appearance of QC due to the bunching 
in a space-continuous structure, on the other hand the 
fermionic system can also describe phenomena of physical 
interest, like the electron diffusion in solid state systems. 

The dynamics of the global system can be described 
by the two-particle Hamiltonian: 

with m indicating the particle mass. Eq. (|18p. that is 
the sum of two Laplacians, is the continuous version of 
the Hamiltonian written in Eq. ([S|) for the CTQW of two 
particles. Both of them are represented, at the initial 
time by a minimum- uncertainty wavepacket, that in 
the space representation, reads: 

0±xo,fe(a^,io) = (^^j expi^ikx J^^j' 

(19) 

where ±xo is the mean position of the space wavepacket 
with variance a, and k = V^^nEk ^[^^ j^^ kinetic energy 
of the particle. Thus the two-fermion (-boson) initial 
wavefunction $/(b) (a;i, 2:2) can be written as: 



(20) 

I 

up of the non-classical correlations. Here Ep is evalu- 
ated by means of the linear entropy e^, in place of the 
von Neumman entropy, in order to make the numerical 
procedure less demanding, as shown in Ref. [27l . Given 
the relation = 1 — Tr(p^), Ep reads 

Ep = Pi(l - Tr(pi)) = Pi (^1 - J dxp\{x,x)J , 

(21) 

where the space-continuous nature of the model is 
explicitly taken into account. In the above ex- 
pression, the probability Pi of finding one parti- 
cle in the Alice side and the other in the Bob 
side, is given by J^^dxi J^°° dx2\^b(f)ixi, X2,t)\'^ , 
while the single-particle density matrix describing 
the subsystem controlled by Alice is pa{x,x') = 
(1/Pi) /+°°da;"$fc(/)(x,x",i)$*(^)(a;',a;",t) with x,x' e 
[-00,0]. 

We examine the model for two different initial condi- 



I 

^f(b)ixi,X2,to) = (/>So,fei(2:i>^o)'/'-So,fe2(2;2,io) - ( + )'feo,fci(2;2,io)'/'-So,fc2(3;i,io), 
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tions: particles with the same velocity, that is fci=fc2, and 
particles with opposite velocities, that is ki—-k2. The 
latter case mimics a collision event between two non in- 
teracting identical particles. 

The former case can be thought as the space- 
continuous analogous of the CTQW on a site lattice with 
two particles initially localized in specific sites. At the 
initial time, both Alice and Bob have a particle: the cen- 
ters of the two single-particle wavefunctions are located 
in different spatial subdomains and the space overlap be- 
tween the wavefunctions is zero. This makes the degree of 
non-classical correlation vanishing, as shown in left panel 
of Fig. [5l As time increases, single-particle wavepackets 
spread out making the spatial overlap between them not 
negligible anymore. Now the probability amplitudes of 
the two-fermion (-boson) wavefunction ^m-){xi,X2) can 
interfere and quantum correlation builds up. The lat- 
ter increases with time ad finally reaches a stationary 
value depending upon the quantum statistics of the par- 
ticles involved: it is higher for bosons. Apart from the 
absence of oscillations which can be related to the space- 
continuous structure, the time evolution of Ep of the 
fermionic (bosonic) system is, at least at short times, in 
qualitative agreement with the dynamics of the quantum 
correlations appearing in the model investigated in the 
previous section. At longer times, due to the different 
boundary conditions of the two systems, a discrepancy 
between the two entanglement-time behavior is observed. 
While in the ring lattice the two-particle wavefunction is 
transmitted through periodic boundary conditions and 
this results in a further increase of the amount of quan- 
tum correlations, the space-continuous structure is sup- 
posed to be infinite and the wavepackets can only spread 
with no additional interference effect due to transmitted 
components. 

For the case of two fermions (bosons) running in op- 
posite directions, QC exhibit a peculiar behavior (see 
the right panel of Fig. [5]). Ep increases while the two 
fermions (bosons) are approaching each other. Specif- 
ically, when the centers of the two wave packets reach 
the minimum distance, two-particle interference gives 
the maximum amount of quantum correlations. Finally, 
Ep drops again to zero once the particles get far apart 
and the corresponding Gaussian wave packets exhibit 
a negligible spatial overlap. This behavior is in agree- 
ment with previous analyses, adopting other entangle- 
ment criteria, of the dynamics of the QC stemming from 
carrier-carrier scattering events in semiconductor nanos- 
tructures [Mlli]. 



V. CONCLUSIONS 

Recent experiments addressed the appearance of non- 
classical correlations in CTWQ of photon pairs in coupled 
waveguide lattices [1, Q ■ Such correlations have been re- 
lated to non trivial interference effects due to the quan- 
tum statistics of the particles. In this perspective, an 



accurate analysis of the quantum correlation created in 
few-particle CTQW models is undoubtedly of great inter- 
est given their experimental feasibility and potential ap- 
plication in quantum communication and quantum com- 
putation theory d, Q . 

Here, we investigated the appearance of QC in CTWQ 
of two non-interacting bosons (fermions) on an ID ring 
lattice with periodic boundary conditions. Specifically, in 
our model the topology of the graphs examined is very 
simple, that is each node is connected to its first neigh- 
bours and this permits to express the Hamiltonian, rul- 
ing the two-particle dynamics, in terms of the so-called 
discrete single-particle Laplacian. Indeed, our CTQW 
well describes the free propagation of two non-interacting 
identical particles in a periodic system. Given the key 
role played by the exchange symmetry into the emergence 
of non-classical correlations, their quantitative evaluation 
in our physical system required the use of a suitable crite- 
rion which takes into account the indistinguishability of 
the particles. To this purpose we adopted the Wiseman 
and Vaccaro approach ^T&\ which allowed estimating the 
degree of quantum correlation between two parties of the 
system (each possessing one particle and accessing to a 
given set of nodes) with no violation of the particle local 
number superselection rules, as occurring in the prepara- 
tion, manipulation and measurement of the experimental 
implementations of CTQW. 

In agreement with the theoretical predictions @, Q, 
results indicate that the building up of QC in our sys- 
tem is due to the two-particle interference effects between 
the propagating wavepackets. Indeed, the time evolu- 
tion of the quantum system affects the degree of corre- 
lation. Specifically, we find that, for CTQW on discrete 
ring lattices with a small number of nodes, the transi- 
tion amplitudes depend sinusoidally on time due to the 
transmission of the wavepackets through periodic bound- 
aries conditions. As a consequence, also QC exhibit a 
periodic-time behavior. They are also affected by the 
quantum statistics of the particles involved in the pro- 
cess and by the partition of the system. Specifically, the 
production of the maximum degree of non-classical cor- 
relation between subsystem accessing non adjacent nodes 
occurs, cyclically, only for two-fermion systems initially 
prepared in a suitable input state. Such a result seems 
consistent with experiments in electron Hanbury Brown- 
Twiss interferometers showing the emergence of QC be- 
tween two identical charge carriers at couples of non ad- 
jacent drains [H [13. For two-particle CTQW on ID 
discrete ring lattice with a large number of nodes, the 
time evolution of the quantum system is less affected by 
transmission through periodic boundary conditions and 
this leads to the disappearance of the cyclical behavior of 
the QC. In particular, the latter are still related to the in- 
terference stemming from the overlap of the wavepackets, 
but they turn out to be almost insensitive to the quantum 
statistics of the particles, a part from small fluctuations 
related to the discrete nature of the ring lattice. 

Finally, we analyzed the building up of QC in the free- 
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propagation of two identical particles in an ID structure. 
Such a system can be considered the continuous space 
limit of an ID ring lattice with a large number of sites and 
vanishing intersite distance. It provides an useful guide- 
line both to validate the results found for the discrete case 
and to compare our outcomes with the analyses of non- 
classical correlations appearing in scattering events in ID 
semiconductor structures ^28, 29]. Two different initial 
conditions have been examined: particles moving with 
the same velocity in the same or in the opposite direc- 
tion. In the former case, we find that the QC of the sys- 
tem first show an initial increase due to the interference 
of the spreading single-particle wavepackets and then, for 
sufficiently long times, reach a stationary value, different 
from zero, depending upon quantum statistics. Such be- 
havior is in qualitative agreement with the one found 
for CTQW in ID ring lattice (apart from the absence 
of the elfects related to the transmission through peri- 
odic boundary conditions and to the discreetness of the 
lattice). On the other hand, for the case of particles prop- 
agating one against the other, as the fermions (bosons) 
get closer, the degree of quantum correlation increases 



and attains its maximum values in correspondence of the 
maximum spatial overlap of the single-particle wavepack- 
ets. When the two particles get away from each other, 
the amount of QC decreases until it vanishes. Such a 
behavior appears to be consistent with the time evolu- 
tion of the entanglement created in binary collisions of 
electrons in semiconductor structures psl l29j . 

Despite the simplicity of the graphs here investigated, 
it seems reasonable to assume that the results obtained 
are representative, at least qualitatively, of the behav- 
ior of the QC between two indistinguishable particles in 
more complex structures, such as 2D lattices or graph 
with larger connectivity. Thus, the analysis here reported 
represents a valuable starting point to investigate the 
emergence of non-classical correlations in other models 
of or fermionic or bosonic CTWQ, where also particle- 
particle interaction could be considered. Furthermore, it 
would be of interest to insert in such models environmen- 
tal noise such as the one due to lattice disorder. Indeed, 
the latter can be viewed as a decoherence source and, as 
a consequence, should result into peculiar phenomena, 
among which the sudden death or the revival of quan- 
tum correlations. 
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FIG. 3. Left panels: two-particle correlation function F^^ 
of the two-fermion CTQW in a ring lattice with N — 70 
evaluated at four different times: "ft = 0, (a), — 5, (b), 
7i = 10, (c), and 7* = 40 (d). The input state is |f f + 
1)/. At short times, fermions exhibit ballistic evolution with 
spatial antibunching. At the final time, as a consequence 
of the interference stemming from the transmission through 
periodic boundary conditions, the two particle wavefunction is 
scattered along the whole lattice. Right panels: two-particle 
correlation function F^^ of the two-boson CTQW in a ring 
lattice with N — 70 evaluated at four different times: 'yt = 0, 
(a), 7t = 5, (b), 7* = 10, (c), and -yt = 40 (d). The input state 
is + l)(j- At short times, bosons propagate ballistically 

along the lattice showing spatial bunching. After reaching the 
boundary, the transmission through the periodic boundary 
conditions induces the spread of two-particle wavefunction all 
over the lattice. 
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FIG. 4. Probability Pi of finding each particle in a differ- 
ent subsystem and von Neumann entropy e„iv, both for two 
bosons and the for two fermions as a function of time for the 
case of a lattice ring with a number of nodes A'' = 50. 
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FIG. 5. (Color online) Left panel: Ep as a function of the 
time for two fermions (dashed line) and two bosons (solid line) 
moving with the same velocity. Such a condition is equiva- 
lent to the one where the particles have zero initial kinetic 
energy. In order to make the numerical implementation of 
wavefunction dynamics simple, here we considered Ek = 0. 
The particles are initially described by two wavepackets with 
initial variance (T=5nm and xo=2Q nm. Right panel: Ep as 
a function of the time for two fermions (dashed line) and two 
bosons (solid line) moving in opposite directions. The parti- 
cles have the same kinetic energy Et = 10 meV, and, at the 
initial time, are also described by two wavepackets with vari- 
ance (T=5nm and xo=2Q nm. In all the numerical calculations 
we take m=9.1x 10"^^Kg. 



